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ABSTRACT 


The  theory  of  dispersive  Rayleigh  waves  coupled  to 
atmospheric  compressionai  waves  is  derived  for  the  case 
of  a  solid  surface  layer.  Numerical  computation  of  phase 
and  group  velocity  curves  indicate  that  an  additional 
branch  may  be  introduced  to  the  dispersion  curves  as  a 
result  of  air  coupling.  Amplitudes  of  waves  propagated 
according  to  the  various  branches  are  briefly  discussed. 


3. 


INTRODUCTION 

Recent  experimental  and  theoretical  work  has  indicated 
that  coupling  of  surface  waves  to  the  atmosphere  is  not  an 
uncommon  phenomena.  Press,  Crary,  Oliver  and  Katz  ^  identi¬ 
fied  air  coupled  flexural  waves  on  a  floating  ice  sheet. 

(2)  (3) 

Press  and  Ewing  made  use  of  Lamb's  elementary  theory  v  7 

to  show  that  coupling  between  air  waves  and  surface  waves  of 
all  types  is  appreciable  when  the  phase  velocity  of  the  surf¬ 
ace  wave  is  very  close  to  the  speed  of  sound  in  air.  These 
authors  also  developed  the  exact  theory  for  air  coupled  flex¬ 
ural  waves  in  a  floating  ice  sheet  and  presented  numerical 
computation  of  dispersion  and  relative  amplitudes. 

An  experimental  investigation  of  air-coupling  to  Rayleigh 
waves  has  recently  been  completed  with  the  cooperation  of  the 
Field  Research  Laboratory  of  Magnolia  Petroleum  Company.  A 
paper  reporting  on  these  results,  entitled  '’Coupling  of  Ground 
Roll  to  Atmospheric  Compressional  Waves"  by  Frank  Press  and 
Maurice  Ewing  will  be  submitted  for  publication  in  "Geophysics" 
shortly . 

Coupling  of  the  atmosphere  with  hydrodynamic  gravity  waves 
is  now  being  investigated. 
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In  this  paper,  the  theory  for  air  coupled.  Rayleigh  waves 
originating  in  an  impulsive  point  source  in  the  air,  is  out 
lined  for  the  case  of  a  solid  surface  layer.  Numerical  com 
putation  of  dispersion  is  made  for  a  particular  case,  the 
results  being  applicable  for  a  source  in  the  air  or  in  the 
ground . 


THEORY 


The  existence  of  air  coupled  Rayleigh  waves  was  shown 


by  Bateman 


(4) 


on  the  assumption  of  an  homogeneous  earth. 


Lee  investigated  the  propagation  of  Rayleigh  waves  in  two 

i-ii  (5>6) 

solid  layers 

We  now  consider  Rayleigh  'Waves  coupled  to  atmospheric 
compressional  waves  for  the  case  of  a  three  layered  space 
formed  by  air  in  the  half  space  z  .  0  to  z  =  -  00,  an 
"upper1’  solid  layer  (1)  between  the  planes  z  =  0  and  z  =  H 
and  a  "lower"  layer  (1)  between  z  =  H  and  z  =  00.  We  assume 
a  point  source  in  the  air  at  a  distance  d  above  the  plane 
z  =  0.  We  denote  by  <*-, ,  ,  ft,  ,  Qz 

the  densities  and  velocities  of  dilatational  and  shear  waves 


in  the  air  and  the  elastic  media  1  and 


(7) 


Applying  the  well  known  method  of  Lamb  v  ,  we  repres¬ 
ent  the  horizontal  component  of  displacement,  q,  and  the 
vertical  component,  w,  by  the  formulae 
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where  r  and  z  are  two  cylindrical  coordinates.  Simple  harmonic 

_  l  cot 


motion  of  the  form 


has  been  assumed.  The  differ- 
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ential  equations  of  motion  may  be  transformed  to  the  following 


system 
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In  order  to  find  the  functions  and  which  satisfy  equat¬ 


ions  (2),  we  must  specify  the  boundary  conditions.  The  stresses 

P„-x*+zr- 


(3) 

At  =  f1  ^  If  J  (4) 

where  A  and  are  Lamb’s  constants,  satisfy  the  conditions 

(P2*),  =L(f?zz)o?  ( P 2 1)[  =  ^  ^  z^O  (j) 


and 


6M*  =(A  4>  (p^)i=(p^),  at  z~~ H  (p) 

If  we  assume  that  no  voids  are  formed  at  the  interface  between 


the  air  and  the  upper  layer,  we  may  write 

Wj-W0  aP  Z  -  0 


(V 


and,  if  a  welded  connection  between  the  two  solid  layers  is  assumed, 
we  have 

WlsW!  ,  %  =  £!>  at  Z‘H.  (3) 

It  is  also  assumed  that  the  displacements  vanish  at  z  =  ±  oo  . 

(8) 

Following  Lamb  ,  we  introduce  two  conditions  which  must  be  satis¬ 
fied  at  the  horizontal  plane  passing  through  the  source  (z  =  d) , 


7.  H.  Lamb,  "On  the  Propagation  of  Tremors  over  the  Surface  of  an 
Elastic  Solid",  Phil.  Trans.  Roy.  Soc.  Lond.,  Ser.  A.  Vol.  203, 
pp.  1  -  42,  1904- 
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7. 


namely,  the  pressure  must  be  continuous 

>  *  z  -  oi 

atZ  z  =  cL  .  ' 
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f*  at-  >"  at 
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and  the  vertical  component  of  the  velocity  is  discontinuous 
at  the  plane  z  -  d 

3*'-  %£=2y7(/n)  (,0) 

The  fimdamental  solutions  of  equations  (1)  satisfying  all 
above  conditions  must  be  written  in  the  well  known  form 

%' = a  **p(y,  z)Jo  (**) ,  M  >1^1 

%”={0&f>(-ytz)  +  CeJtp(jf,z)J^ (*!t) ,  /ct/i  jzf  s.  o 
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where  JQ  (k£)  is  the  Bessel  function  and  the  positive  real  or 
positive  imaginary  constants  /;  are  given  by  relations 
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in  order  to  satisfy  equations  (l) .  Nine  arbitrary  constants. 


(”) 


A,B,C,...  .  P  figuring  in  the  expressions  (Il)-(l6)  may  be 
chosen  so  that  the  nine  conditions  ( 5) , (6) , (7) , (8) , (9)  and  (10) 
are  satisfied. 

On  substituting  the  functions  (11) -(16)  in  these  boundary 
conditions,  we  obtain  the  following  system  of  equations  to  deter¬ 
mine  the  values  of  A,B,C,....P: 
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Two  of  the  coefficients  A,B,....P  may  be  easily 
eliminated.  From  the  two  last  equations  it  follows 

'/R  ^  o  r  \ 
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Therefore,  the  equations  (18)-(a4)  take  the  form 
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9. 


Since  the  determinant  of  these  seven  equations  is 


not  identically  equal  to  zero  we  can  find  definite  values 
of  the  coefficients  A,B,...P.  We  can  then  apply  the  method 


velocity  using  the  Fourier-Bessel  integral  by  taking  Y~kdk 
and  integrating  with  respect  to  k  from  0  to  00.  Then  we 
can  obtain  solutions  satisfying  all  conditions  including  the 
requirements  for  a  point  source.  It  may  be  difficult  to 
evaluate  the  integrals  representing  the  solutions,  since 

the  quantities  A,B . P  are  complicated,  owing  to  the  high 

order  of  the  determinant. 

We  restrict  ourselves  at  this  point  to  the  analysis 
of  the  frequency  equation  which  is  given  by  the  condition 
that  the  above  determinant  is  equal  to  zero. 

Since  the  coefficients  a,d,....  of  the  system  (19) 
depend  on  the  elastic  constants  and  on  the  quantities  pL  } 


(Xc  }  jo;  and  £  which  may  be  expressed  as  functions 

of  the  phase  velocity,  C  -  00 /k.  ,  this  determinant  may 


be  considered  as  a  function  of  the  phase  velocity  c  and 


of  the  frequency 


.  Equating  the  determinant  to 


zero,  we,  therefore,  obtain  an  equation  which  defines  an 
implicit  relationship  between  frequency  and  the  phase  veloc¬ 
ity. 


9.  H.  Lamb,  "On  Propagation  of  Tremors 
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Rewriting  the  frequency  equation,  we  omit  all  factors 
which  can -not  be  equal  to  zero  and  transform  the  remaining 
determinant  to  reduce  the  number  of  elements.  After  these 
transformations,  the  frequency  equation  takes  the  form: 
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Developing  the  determinant,  we  put 

c 


A 


and 


n 


n. 


Introducing  Lee’s  notations 

X  -  f  v‘-  -  IJ 


10 


y  =  yV^-/J  =  ^V 


w 


Z  =^^-^-2(^-0  =  X-\ 


we  put 


<j, -xz-v.w ,  yv,  £  = 
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The  frequency  equation  then  takes  the  form 
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10.  A.  W.  Lee,  “The  Effects  of  Geologic  Structure . “,  and 

“Further  Effects  . . " 


having  the  factors 


The  last  terms  of  the  coefficients 


l 


and  /nc  represent  the  influence  of  the  air. 

For  numerical  computations,  it  was  assumed  that  Poisson's  ratio 
£3  =0.25  for  both  solid  layers,  i.e.  and  f^Z  4 

The  following  data  were  used: 

<xc=lo7o  ■ft/szc  ,  ji,  =  8oo 


Si  =  /.J?  ,  SS  ~  O.OC!  ,  ^  =  /3.77 . 

f>  f,  ’ 

These  values  were  chosen  as  most  representative  of  near  surface 
conditions  often  encountered  in  seismic  prospecting.  It  is 
obvious  that  air  coupling  effects  are  negligible  when  ,  >°^0 

since  c  >  for  this  case,  and  the  air  coupling  term  of 
the  frequency  equation  which  contains  the  factor  <?0 / n0 
is  very  small. 

In  the  case  of  air  coupled  Rayleigh  waves,  real  values  of  kH 
must  correspond  to  the  interval 


.  9/97  6^  =  t 

A 

The  upper  limit  is  equal  in  this  case,  to  1.3375  (-  )• 

For  larger  values  of  V,  the  roots  of  the  frequency  equation 
become  complex.  We  have  computed  the  real  and  imaginary  parts 
of  these  roots.  They  are  given  in  Table  II.  All  real  roots 
computed,  and  the  real  parts  of  complex  roots,  are  given  in  Table 
I.  The  roots  were  computed  by  successive  approximations,  using 
Lee's  computations  to  obtain  the  first  approximation.  It  may  be 
remarked  that  the  values  of  L/H  as  given  by  Lee  to  two  or  three 
figures,  lead  to  values  of  kH  which  may  differ,  by  some  units  of 


the  third  figure,  from  the  more  exact  values  computed  with 
five  figures  directly  from  the  frequency  equation.  Group 
velocity  U  was  obtained  by  graphical  differentiation  of  the 


phase  velocity  curve,  using  the  familiar  equation 

■ 

The  phase  velocity  curve  in  the  case  of  air  coupled 
Rayleigh  waves  differs  very  little  from  the  curve  given  by 
Lee  until  the  neighbourhood  of  the  critical  point  V  =  1.3375 
is  reached  (Figure  I) .  At  this  point,  it  deviates  to  the 
left  and  intersects  the  V-axis  at  a  point 


The  real  part  of  complex  roots  above  the  critical  value 
is  very  close  to  the  corresponding  roots  of  the  frequency  equa¬ 
tion  without  air  connected  terms.  But  it  deviates  to  the 
right  in  the  neighbourhood  of  the  critical  value  and  "ends" 
at  the  point  a7//  -  ^/.  25~/5~  of  the  line  V  -  1.3375* 

Results  of  these  computations  are  presented  in  dimen¬ 
sionless  form  in  Figure  I  where  and  are  plotted 

against  kH.  The  heavy  lines  indicate  where  atmospheric  influ¬ 
ence  is  negligible  and  the  results  are  similar  to  those  obtained 
by  Lee.  The  dashed  lines  represent  new  branches  introduced  by 
air  coupling. 

It  is  well  known  that  the  velocities  and  periods  of  waves 
propagated  through  a  dispersive  medium  can  be  obtained  from 
the  corresponding  group  velocity  curve.  To  obtain  the  ampli¬ 
tudes  of  waves  corresponding  to  the  different  branches  of  the 
group  velocity  curve,  one  must  compute  excitation  functions  which 


depend  upon  the  nature  and  position  of  the  source,  the  location 
of  the  receivers,  and  implicitly  on  the  roots  of  the  frequency 
equation.  Although  the  mathematical  procedure  is  straight¬ 
forward,  the  numerical  work  required  is  extremely  lengthy.  We 
can,  therefore,  discuss  the  frequencies  and  velocities  of 
Rayleigh  waves  coupled  to  the  atmosphere,  but  the  relative 
amplitudes  of  waves  for  different  source  and  receiver  locations 
cannot  be  given.  Existing  calculations  for  air  coupling  to 
another  type  of  surface  wave  can  be  used,  however,  to  infer 
something  about  relative  amplitudes. 

The  group  velocity  curve  in  Figure  I  is  divided  into 
branches  I,  II,  III,  each  of  which  represents  a  different 
train  of  waves.  Branch  I  corresponds  to  the  dispersive 
train  of  Rayleigh  waves  observed  on  seismograms  of  earthquakes. 
Branch  I  also  accounts  for  the  dispersive  Rayleigh  waves  usu¬ 
ally  associated  with  ground  roll.  These  waves  first  appear  as 
long  period  arrivals  travelling  with  the  speed  of  Rayleigh 
waves  in  the  bottom  layer,  U  -  .9194^  •  Succeeding  waves 

gradually  decrease  in  period  since  kH  increases  as  the  group 
velocity  decreases.  Waves  continue  to  arrive  with  decreasing 
period  until  a  time  corresponding  to  propagation  at  the  min¬ 
imum  value  of  group  velocity.  From  the  analogous  theory  re¬ 
ferred  to  above,  we  would  expect  waves  propagated  according 
to  branch  I  to  predominate  for  a  source  located  within  the 
solid  layers  recorded  by, pickup  on  the  solid  layer.  Waves 
with  group  velocity  values  near  the  minimum  in  Branch  I  have 


phase  velocities  approaching  the  speed  of  sound  in  air 
from  the  side  c  >  cxo  .  These  waves  are  attenuated  since 
kH  is  complex  and  has  an  increasingly  large  imaginary  com¬ 
ponent  as  c  ->  oc0  (Table  II) . 

No  significance  should  be  attached  to  the  occurrence  of  a 
minimum  value  of  group  velocity  near  the  critical  point 
c  —  VQ  .  Had  we  chosen  other  values  for  the  elastic  con¬ 
stants  of  the  layers,  the  phase  velocity  corresponding  to 
the  minimal  group  velocity  could  have  been  much  different 
from  VQ  . 

Branch  II  represents  a  dispersive  train  of  waves  be¬ 
ginning  as  a  high  frequency  arrival  at  a  time  corresponding 
to  propagation  with  the  velocity  of  Rayleigh  waves  in  the 
surface  layer.  The  frequency  of  these  waves  decreases  as 
time  progresses  until  the  time  corresponding  to  propaga¬ 
tion  at  the  minimum  group  velocity  when  the  waves  of  Branch 
I  and  Branch  II  merge  to  form  a  single  train  of  waves  hav¬ 
ing  a  discrete  frequency. 

Branch  III  represents  an  additional  train  introduced 
by  coupling  of  Rayleigh  waves  to  atmospheric  compressional 
waves.  This  train  begins  at  a  time  corresponding  to  propa¬ 
gation  at  the  speed  of  sound  in  air  and  continues  with  al¬ 
most  constant  frequency  until  the  time  t  =  r/.44£j_  *  The 
phase  velocity  of  these  waves  should  be  close  to  the  speed 
of  sound  in  air.  From  experimental  and  theoretical  results 
referred  to  earlier,  we  might  expect  these  waves  to  be  prom- 


inent  for  a  source  in  the  air  recorded  by  a  pickup  on  the 
ground,  and  for  a  source  within  the  ground  recorded  by  a 
microphone  in  the  air.  There  is  an  additional  branch  with 
dispersion  features  similar  to  Branch  III,  corresponding 
to  the  complex  phase  velocities  whose  real  parts  are  greater 
than,  but  close  to  oc^ 

Higher  modes  of  propagation  exist  at  correspondingly 
higher  frequencies. 


COHCLUSIQHS 

The  existence  of  a  surface  layer  in  which  the  speed 
of  Rayleigh  waves  is  less  than  the  speed  of  sound  in  air 
introduces  a  new  branch  to  the  dispersion  curves,  correspond¬ 
ing  to  a  train  of  approximately  constant  frequency  waves  be¬ 
ginning  with  the  arrival  of  the  air  wave  and  continuing  to  a 
time  depending  on  the  distance  and  the  elastic  constants  of 
the  layers.  It  is  probable  that  these  new  waves  are  espec¬ 
ially  prominent  for  a  source  within  the  air  recorded  with  a 
seismograph  and  for  a  source  within  the  ground  recorded  with 
a  microphone. 
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1.  This  is  a  more  exact  value  than  that  obtained 
from  Lee's  data  (1.086) 
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Real  part 
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U  Figure  1.  Phase  and  group  velocity  curves  for  air  coupled  Rayleigh  waves. 
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